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Abstract 

The  goal  of  this  study  is  to  develop  a  system-level  dynamic  model  of  a  proton-exchange  membrane  (PEM)  fuel  cell  that  is  capable  of 
characterizing  the  mixed  effects  of  temperature,  gas  flow,  and  capacitance,  with  particular  emphasis  focused  on  system  transient  behavior. 
The  fuel  cell  system  is  divided  into  three  control  volumes  and  thus  a  lumped-parameter  model  for  these  sub-systems  is  established  using  a 
combination  of  intrinsic  mechanistic  relations  and  empirical  modeling.  The  dynamic  model  is  simulated  using  SIMULINK.  The  analysis 
illustrates  the  complicated  dynamic  interactions  between  various  components  and  effects  within  a  fuel  cell  system,  and  demonstrates  the 
necessity  of  the  proposed  approach  of  separate  control  volumes.  Numerical  studies  are  correlated  to  a  single-cell  experimental  investigation, 
and  a  protocol  for  parameter  identification  is  explored  to  refine  the  model  fidelity.  The  proposed  fuel  cell  model  can  accurately  predict  the 
dynamic  behavior  and  exhibits  excellent  agreement  with  experimental  results.  This  model  can  be  readily  employed  in  the  optimization  and 
real-time  control  of  PEM  fuel  cells  installed  in  practical  automotive  or  stationary  applications. 
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1.  Introduction 

A  proton-exchange  membrane  fuel  cell  (PEMFC)  is  an 
electrochemical  device  that  converts  the  chemical  energy 
of  hydrogen  and  oxygen,  with  the  aid  of  electro-catalysts, 
directly  into  electrical  energy.  After  four  decades  of  re¬ 
search  and  development,  this  device  has  reached  the  test  and 
demonstration  phase  [1].  Typically,  the  analysis  and  design 
of  PEMFC s  are  centered  around  the  membrane-electrode  as¬ 
sembly  (ME  A).  The  modeling  of  an  ME  A  involves  the  char¬ 
acterization  of  physical  environment  of  the  electrochemical 
reaction,  the  transport  phenomena  of  gas  (hydrogen,  oxy¬ 
gen,  water  vapor,  etc.),  water,  proton  and  current,  and  the 
relationship  between  the  fuel  cell  voltage,  current,  temper¬ 
ature,  material  (electrode,  catalyst  and  membrane)  proper¬ 
ties  and  transport  parameters.  Traditional  ME  A  modeling 
is  mechanistic  in  essence  [2-7].  Attempts  have  been  made 
to  investigate  multi-species  diffusion  through  the  substrate 
and  the  diffusion  layer  of  the  electrode,  the  reaction  kinetics 
in  the  active  layer  (catalyst),  and  proton  and  water  trans¬ 
port  through  membrane.  The  resultant  models  are  generally 
governed  by  a  set  of  complex  partial  differential  equations. 
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While  the  mechanistic  models  are  built  upon  the  rigorous 
mathematical  descriptions  of  fundamental  physics,  they  nor¬ 
mally  have  very  complicated  expressions  and  oftentimes  cer¬ 
tain  key  physical  parameters  are  difficult  to  quantify,  which 
may  in  turn  reduce  the  modeling  effectiveness.  Empirical 
modeling,  by  mapping  the  fuel  cell  voltage  as  a  function 
of  various  contributing  variables  [8-11],  on  the  other  hand, 
has  certain  advantages  in  some  applications.  Although  this 
approach  does  not  usually  examine  in  depth  the  underly¬ 
ing  electrochemical  process,  it  utilizes  the  phenomenologi¬ 
cal  description  based  on  benchmark  experimental  study  to 
reduce  the  model  complexity  as  well  as  computational  time. 
Typically,  such  empirical  MEA  sub-models  are  used  in  large 
fuel  cell  system-level  models  [12-15]. 

While  MEA  modeling  constitutes  the  kernel  of  PEMFC 
study,  the  dynamic  modeling  and  simulation  at  system  level 
are  even  more  complicated  and  involve  fluid-structure-heat- 
electrochemistry  coupled  behavior.  Amphlett  et  al.  [12]  pro¬ 
posed  a  fuel  cell  dynamic  model  which  includes  fuel  cell 
body  temperature  variation  effect  but  used  steady-state  flow 
rate  governing  equations.  Pukrushpan  et  al.  [13]  developed 
a  system-level  model  that  includes  compressor,  supply  and 
return  manifolds,  humidifier,  and  anode  and  cathode  chan¬ 
nels.  Similar  work  was  performed  by  Iqbal  [14]  who  stud¬ 
ied  a  hybrid  energy  system  with  an  embedded  fuel  cell. 
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Yerramalla  et  al.  [15]  considered  the  humidifier  and  stack 
pressure  in  the  system  level  fuel  cell  modeling.  It  is  worth 
noting  that  while  the  aforementioned  studies  have  laid  down 
a  foundation  for  fuel  cell  system-level  modeling,  they  had 
their  respective  emphases  on  certain  aspects  and  have  not 
addressed  the  complete  dynamic  effects  within  the  system. 
A  common  simplification  in  references  [13-15]  is  that  the 
temperature  transient  behavior  was  neglected.  On  the  other 
hand,  it  has  been  well  recognized  that  temperature  is  a  key 
parameter  affecting  the  performance  of  the  fuel  cell.  While 
low  temperature  of  the  cell  could  lead  to  flooding  of  the 
membrane,  high  temperature  could  cause  a  dry  membrane, 
both  of  which  could  reduce  the  performance  of  the  fuel  cell 
in  terms  of  heat  and  water  management.  Indeed,  the  fuel  tem¬ 
perature  and  gas  flow  play  key  roles  in  fuel  cell  operation, 
and  their  dynamic  effects  should  be  simultaneously  and  ac¬ 
curately  characterized  in  a  fuel  cell  model  in  order  for  such 
a  model  to  be  effective  within  the  entire  operating  range.  It 
should  be  noted  that  some  of  the  previous  studies  have  not 
been  completely  correlated  to  experimental  investigations. 

It  is  the  goal  of  this  study  to  develop  a  system-level  dy¬ 
namic  model  of  a  PEMFC  that  incorporates  the  complicated 
temperature,  gas  flow  and  capacitance  effects  under  operat¬ 
ing  conditions.  The  focus  is  on  the  dynamic  and  transient 
properties  of  the  system.  We  will  use  the  control  volume 
approach  [16]  to  develop  a  set  of  complete  dynamic  equa¬ 
tions  that  govern  the  system  dynamics.  The  fuel  cell  system 
is  divided  into  three  control  volumes,  the  anode  channel,  the 
cathode  channel,  and  the  fuel  cell  body.  For  each  control 
volume,  the  establishment  of  a  lumped-parameter  dynamic 
model  is  realized  using  a  combination  of  intrinsic  mech¬ 
anistic  relations  and  empirical  modeling.  The  system-level 
model  is  simulated  using  SIMUFINK,  and  a  series  of  anal¬ 
yses  are  carried  out.  Numerical  studies  are  then  correlated 
to  single-cell  experimental  investigations,  and  a  protocol  for 
parameter  identification  is  explored  to  refine  the  model  fi¬ 
delity.  This  model  is  capable  of  characterizing  the  transient 
dynamic  properties  of  a  fuel  cell  system,  and  can  be  read¬ 
ily  employed  in  the  optimization  and  real-time  control  of 
PEMFC s  installed  in  practical  automotive  or  stationary  ap¬ 
plications. 


2.  Fuel  cell  modeling  using  control  volume  approach 

In  this  study,  we  develop  a  mathematical  model  to  investi¬ 
gate  the  mass  and  energy  dynamic  interactions  as  well  as  the 
electrochemical  reaction  within  a  PEMFC  by  using  lumped 
parameters.  As  shown  in  Fig.  1,  hydrogen  enters  the  anode 
channel,  part  of  which  diffuses  through  the  anode  diffusion 
layer  to  the  active  catalyst  layer,  where  the  hydrogen  disso¬ 
ciates  into  protons  and  electrons.  The  protons  travel  from  the 
anode  catalyst,  through  the  membrane,  to  the  cathode  cat¬ 
alyst.  In  the  meantime,  oxygen  enters  the  cathode  channel, 
part  of  which  diffuse  through  the  cathode  diffusion  layer  to 
the  active  catalyst,  where  the  oxygen  dissociates  and  com- 
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Fig.  1.  Schematic  illustration  of  a  PEM  fuel  cell. 


bines  with  protons  and  electrons  to  form  water,  while  heat 
and  current  are  generated.  The  heat  produced  during  the  re¬ 
duction  reaction  process  will  warm  up  the  fuel  cell  body, 
which  will  affect  the  mass  diffusion  within  the  anode  dif¬ 
fusion  layer,  cathode  diffusion  layer  and  membrane.  Typi¬ 
cally,  the  heat  transfers  to  ambient  and  gas  flow  channels 
through  convection  and  radiation.  The  resulting  gas  temper¬ 
ature  variation  will  change  the  pressure  and  flow  rate  within 
the  channels,  which  therefore  affects  the  gas  diffusion  within 
the  diffusion  layers.  All  these  dynamic  interactions  will  have 
a  direct  effect  on  the  fuel  cell  output  voltage. 

Throughout  this  paper,  we  assume  an  ideal  gas  mixture 
[19],  which  means  that  each  mixture  component  behaves  as 
an  ideal  gas  as  if  it  were  alone.  The  momentum  equation  is 
not  included,  as  we  are  concerned  with  macroscopic  level 
dynamic  properties.  A  no- slip  condition  at  the  fixed  wall  is 
assumed,  and  the  viscous  normal  stresses  at  the  inlet  and 
outlet  are  neglected. 


2.7.  Modeling  of  mass  and  energy  transfer  using  control 
volume  approach 


2.1.1.  General  governing  equations 

In  this  subsection,  we  present  the  basic  governing  equa¬ 
tions  for  mass  and  energy  interactions.  The  detailed  descrip¬ 
tion  of  the  mathematical  model  for  each  control  volume  will 
be  given  in  the  subsequent  subsections. 

For  an  arbitrary  control  volume,  the  continuity  equation 
in  integral  form  is  given  as  [16] 


pdV  + 


# 

c.s. 


p(V  •  n)  d  A  =  0 


and  the  energy  equation  in  integral  form  is  [16] 


d 

dt 


J I  j  pet  dV  +  ph(V  •  n)  d A 

c.v.  c.s. 


d  Q  dW 
dt  dt 


where  p  is  the  material  density  within  the  control  volume, 
V  the  velocity  vector  of  flow  fluid,  h  the  unit  vector  normal 
to  the  control  surface  toward  the  outside,  et  the  specific 
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internal  energy,  h  the  specific  enthalpy,  Q  the  heat  added 
to  the  control  volume,  and  W  the  work  done  on  the  control 
volume.  Note  that  typical  gas  flow  velocity  is  low  under  fuel 
cell  operating  conditions.  The  kinetic  and  potential  energy 
of  the  gases  is  thus  neglected  in  the  energy  equation. 

In  a  fuel  cell  system,  there  are  complicated  mass  and  en¬ 
ergy  interactions  between  the  gas  flow  channels,  cell  body 
and  ME  A.  In  order  to  quantify  such  interactions,  in  this  pa¬ 
per  we  define  three  control  volumes,  the  anode  channel,  the 
cathode  channel,  and  the  body  of  the  fuel  cell.  The  anode 
channel  control  volume  is  the  lumped  volume  from  the  an¬ 
ode  inlet  to  the  anode  outlet.  The  cathode  channel  control 
volume  is  similarly  defined  at  the  cathode  side.  The  MEA 
and  all  other  metal  structures  involved  in  the  fuel  cell  will 
form  the  fuel  cell  body  control  volume.  When  a  fuel  cell 
stack  is  considered,  the  anode  channel,  cathode  channel  and 
body  control  volumes  will  be  the  summation  of  correspond¬ 
ing  control  volumes  of  each  single  cell.  The  physical  prop¬ 
erties  associated  with  the  control  volumes  are  assumed  to  be 
uniform.  That  is,  in  this  study  within  the  control  volume  we 
consider  the  average  effect.  In  the  subsequent  subsections, 
we  will  use  Eqs.  (1)  and  (2)  to  derive  detailed  governing 
equations  corresponding  to  the  three  control  volumes.  Since 
the  characterization  of  water  vapor  is  extremely  complicated 
and  involves  phase  change,  such  effects  are  not  considered 
in  the  present  study. 


2.7.2.  Anode  and  cathode  channels 

For  the  control  volumes  corresponding  to  the  anode  and 
cathode  channels,  the  first  term  in  Eq.  (1)  becomes  the  mass 
change  rate  of  the  respective  gases.  The  second  term  in 
Eq.  (1)  is  the  mass  flow  rate  through  the  surface  of  the 
control  volumes.  For  an  individual  channel,  three  kinds  of 
mass  flow  rate  are  involved,  which  correspond  to,  the  gas 
entering  the  flow  channel  from  the  inlet,  the  gas  leaving  the 
channel  from  the  outlet,  and  the  gas  diffusion  through  the 
anode/cathode  diffusion  layers,  respectively,  as  described  in 
Fig.  1 
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c.s. 


p(V  •  n)  d A 


// 


inlet 


// 


p(V  •  h)  d A  +  p(V  •  n) d A 


// 

electrode 


outlet 


p(V  •  n) d A 


The  velocity  profiles  are  required  for  the  calculation  of 
Eq.  (3).  In  general,  the  instant  velocity  profile  can  be  ob¬ 
tained  by  solving  the  momentum  equation  simultaneously. 
This  paper  deals  with  the  macroscopic  level  dynamic  behav¬ 
ior  of  the  fuel  cell,  and  the  mass  flow  through  the  inlet/outlet 
control  surfaces  will  be  simplified  as  the  traditional  nozzle 
flow  rate  equation  [13] 


p(V  •  n)  d A  =  k(A p) 


where  A p  denotes  the  pressure  difference,  and  k  the  mass 
flow  rate  coefficient  which  in  turn  depends  upon  temperature 
and  pressure  difference. 

The  third  term  on  the  right-hand  side  of  Eq.  (3)  in¬ 
volves  a  complicated  process.  Hydrogen  diffuses  through 
the  porous  gas-diffusion  anode  to  anode  catalyst,  where  the 
electrochemical  reaction  takes  place  to  form  a  proton  and  an 
electron;  oxygen  diffuses  through  the  porous  gas-diffusion 
cathode  to  cathode  catalyst,  where  the  reduction  reaction 
takes  place  to  form  water  and  heat.  During  the  electro¬ 
chemical/reduction  reaction,  current  is  generated.  Although 
such  processes  at  the  catalyst  layer  are  complicated  and  still 
subject  to  intensive  study,  it  complies  with  the  macroscopic 
level  mass  and  species  conservation  law,  i.e.,  the  consumed 
hydrogen  and  oxygen  have  a  definitive  relationship  with  the 
current  generated  [13] 


p(V  -n)dA  =  N  —  M 

nF 


electrode 


where  i  denotes  the  current,  N  the  number  of  cells  in  a 
stack,  and  F,  M  and  n  are  Faraday  constant,  gas  mole  mass 
and  constant  corresponding  to  different  gas  species,  respec¬ 
tively. 

The  general  energy  equation  (Eq.  (2))  can  be  applied  to 
the  anode  and  cathode  channels  in  a  similar  manner.  The 
first  term  on  the  left-hand  side  of  Eq.  (2)  represents  the  in¬ 
ternal  energy  change  rate  within  the  control  volume.  The 
second  term  in  Eq.  (2)  is  the  heat  transfer  rate  due  to  the 
mass  transport  through  the  control  surface.  Such  mass  trans¬ 
port  involves  the  hydrogen  in  the  anode  channel  that  diffuses 
through  the  porous  anode  diffusion  layer  to  the  cell  body 
and  the  oxygen  in  the  cathode  channel  that  diffuses  through 
the  porous  cathode  diffusion  layer  to  the  cell  body.  On  the 
right-hand  side  of  Eq.  (2),  the  first  term  ( dQ/dt )  represents 
the  rate  of  heat  added  to  the  control  volume,  which  is  gener¬ 
ally  caused  by  convection  heat  transfer  or  possible  radiation 
between  the  cell  body  and  the  channels  and  between  the  cell 
body  and  the  surrounding  ambient.  As  the  temperature  of 
a  PEMFC  body  is  relatively  moderate,  the  radiation  effect 
can  neglected.  Newton’s  cooling  law  governs  the  convection 
heat  transfer  mechanism 


d Q 

d  t 


—  /ZrnnvA 


conv^conv 


AT 


where  hcom  is  the  convection  heat  transfer  coefficient,  Aconv 
the  convection  surface  area,  and  AT  denotes  the  temperature 
difference  between  the  cell  body  and  the  channels/ambient. 
In  this  paper,  the  viscous  normal  stresses  at  the  inlet,  outlet 
and  electrode  interface  of  mass  diffusion  are  neglected,  and 
the  work  done  on  the  control  volume  ( dW/dt )  disappears. 


2.1.3.  Cell  body 

It  is  well  known  that  most  parts  of  the  cell  body  are  com¬ 
posed  of  metal  structures.  The  continuity  equation  (Eq.  (1)) 
for  the  cell  body  control  volume  is  only  concerned  with 
the  species’  conservation  of  electrochemical/reduction  re- 
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action.  That  is,  water  is  formed  through  the  electrochemi¬ 
cal/reduction  reaction  by  consuming  hydrogen  and  oxygen 

Ph2(E  •  n)dA 

electrode,  an 


+ 


// 


Po2(V  -n)dA  =  mn2o 


(7) 


electrode,  ca 


Hereafter,  the  variable  m  represents  the  mass,  and  the  sub¬ 
script  designates  specific  species. 

The  energy  equation  (Eq.  (2))  for  the  cell  body  control 
volume  is  now  considered.  The  first  term  on  the  left-hand 
side  is  the  cell  body  internal  energy  accumulation  rate.  The 
second  term  is  the  heat  transfer  to  the  cell  body  via  mass 
transport.  As  discussed  in  the  Section  2.1.2,  such  mass  trans¬ 
port  includes  the  diffusion  of  gases  through  the  electrodes 
and  the  water  produced  by  the  electrochemical/reduction  re¬ 
action  that  flows  out  from  the  cell  body.  There  are  several 
heat  transfer  effects  involved  in  the  term  dQ/dt ,  which  in¬ 
clude  the  heat  added  to  the  cell  body  through  the  electro¬ 
chemical/reduction  reaction  as  well  as  the  heat  transferred 
to  the  ambient  and  the  channels  via  convection 


d Q 

dt 


—  Gconv,an  T  2conv,ca  T  <2conv,amb 


+  A  Hr 


■T  jj  PUJV 


n)  d A 


eletrode 


where,  specifically,  <2Conv,an  is  the  heat  convection  between 
the  cell  body  and  the  anode  channel,  <2Conv,ca  the  heat  con¬ 
vection  between  the  cell  body  and  the  cathode  channel, 
Gconv,amb  the  heat  convection  between  the  cell  body  and  the 
ambient,  and  A//r  t  the  lower  heating  value  of  the  fuel.  All 
the  heat  convection  terms  here  are  assumed  to  comply  with 
Newton’s  cooling  law  as  shown  in  Eq.  (6).For  the  cell  body, 
the  term  dW/dt  represents  the  electric  power  converted  from 
electrochemical/reduction  reaction,  i.e. 


dW 

dt 


—NV celP 


where  Vceii  denotes  the  cell  voltage  that  will  be  discussed  in 
detail  in  Section  2.2. 


2.1.4.  Dynamic  equations  of  three  control  volumes 

Based  upon  the  above  discussions,  we  derive  detailed 
lumped-parameter  dynamic  equations  governing  all  three 
control  volumes.  Here  we  assume  uniform  thermo-physical 
properties  throughout  the  procedure. 

For  the  anode  channel,  the  hydrogen  mass  change  rate  is 
such  that 


ilflm>dv= 

c.v. 


dmn2 

dt 


(10a) 


and  the  hydrogen  mass  transport  rate  through  the  anode 
channel  control  surface  is  given  as 


-O  pu2(V  -n)dA 

c.s. 

=  ^up,an  (Ps  ,an  Pan)  ^down,an  (P  an  Patm) 


—  N  —  Mh? 
IF  2 


(10b) 


where  &up,an  and  £up?an  are  inlet  and  outlet  mass  flow  rate 
coefficients,  and  ps,an,  Pan  and  pa tm  represent,  respectively, 
the  hydrogen  source  pressure,  the  pressure  within  the  anode 
control  volume  and  ambient  pressure.  The  internal  energy 
change  rate  of  hydrogen  within  the  anode  channel  control 
volume  is  given  as 


Pft2etdV  = 


d(mu2cv,ii2Tm) 

dt 


(10c) 


where  cv^2  Is  the  specific  heat  of  hydrogen  with  constant 
volume,  and  Tan  the  hydrogen  temperature  within  the  an¬ 
ode  channel  control  volume.  The  heat  transfer  to  the  anode 
channel  control  volume  through  the  hydrogen  mass  transfer 
across  the  control  surface,  which  includes  inlet,  outlet  and 
electrode  diffusion  layer,  is  given  as 


(V  •  n)  dA 


up,  an  (Ps  ,an 
"  ^down,an  (P 


an 


Pan)C/?,H2  (7 an  ,in 
—  Patm)Cp,H2  (^an 


-  N  —  Mu2cp^2(Tm 


(10d) 


where  Pan,in  and  T0  represent  the  flow-in  hydrogen  temper¬ 
ature  and  the  reference  temperature,  respectively.  The  con¬ 
vection  heat  transfer  from  the  cell  body  to  the  anode  channel 
is  described  as 


d Q 

dt 


—  ^conv,  anybody  ^an) 


(10e) 


where  &conv,an  —  ^conv,an  AConv,an->  ^conv,an  IS  the  heat  COn 
vection  coefficient,  and  Aconv,an  the  surface  area  of  heat  con¬ 
vection. 

We  can  derive  similar  relations  for  the  cathode  chan¬ 
nel.  The  following  equations  are  parallel  to  those  given  in 
Eqs.  (10a)-(10e),  and  the  relevant  subscripts  are  replaced  by 
“ca”  (cathode)  and  “O2”  (oxygen) 


llfh’-iv= 

C.V. 


dmo2 

dt 


(Ha) 


Po2(V  ■  n)  dA 

c.s. 

=  ^up,ca(Ps,ca  —  Pea)  —  ^down,ca(Pca  —  Patm) 


d 

dt 


-N  —  Mo, 

4  F  2 

///  p°2  e>  d  l/  = 
c.v. 


d(mo2cV'02  7ca) 
dt 


(lib) 

(11c) 
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-  O  po2ho2(V  ■  n)dA 

J  J 

C.S. 

=  ^up,ca(Ps,ca  —  Pc&)cp,02  (^ca,in  —  ^o) 

—  ^down,ca(Pca  —  Patm)Cp,02  (^ca.  —  ^o) 

N  Mo2c  po2  (Tca  T0) 

(lid) 

dQ 

^  —  ^conv,ca(7body  Fca) 

(lie) 

Finally  we  consider  the  cell  body.  Since  most  parts  of  the 
body  are  composed  of  metal  structures  and  MEA  and  gas 
diffusion  and  water  produced  have  little  effect  on  the  total 
mass  change  of  cell  body,  we  thus  have 


d 

Jt 


c.v. 


Pbody  dF  —  0 


(12a) 


The  mass  transport  rate  through  the  cell  body  control  surface 
(MEA)  is  given  as 


Pbody ( V  -n)dA  =  N— Mh2  +  N—Mq2  -  mu2o 


c.s. 


(12b) 


Note  the  above  equation  includes  the  gas  diffusion 
through  the  electrodes  and  water  formed  by  electrochemi¬ 
cal/reduction  reactions.  The  internal  energy  change  rate  of 
cell  body  is  such  that 


d 

Jt 


c.v. 


,  dTbody 

PbodyGdv  — ^bodyCp,body  ^ 


(12c) 


where  embody  is  the  average  specific  heat  of  the  cell  body. 
The  heat  added  to  and  the  work  done  on  the  cell  body  control 
volume  can  be  characterized  as 


d Q 

d  t 


+ 


d  W 
d  t 


^conv,an  (2an  ^body)  T  ^conv,ca  (Tca 
T  ^conv,amb  (  Aimb  —  Aiody) 


+  AHRjN  —  Mu2  -NV^i 
Zr 


^body) 

(12d) 


Such  effects  include  the  convection  heat  transfer  from  the 
surrounding  structures  (the  first  three  terms  at  right-hand 
side),  heat  produced  by  the  electrochemical/reduction  reac¬ 
tion  (the  fourth  term  on  the  right-hand  side),  and  converted 
electric  power  (the  last  term). 


where  Pu2  is  the  hydrogen  partial  pressure,  Po2  the  oxygen 
partial  pressure,  Ts  the  body  temperature,  i  the  load  current, 
and  a  represents  other  possibly  related  parameters.  In  this 
study,  we  use  the  same  approach  to  perform  the  electrochem¬ 
ical  modeling,  which  is  based  upon  a  curve-fitting  method 
[9] .  It  is  worth  noting  that,  in  classical  empirical  MEA  mod¬ 
eling  [9-11,17],  the  effective  surface  partial  pressures  are 
used  in  Eq.  (13).  In  this  study,  in  order  to  simplify  the  em¬ 
pirical  MEA  modeling  for  large  fuel  cell  system,  we  use  the 
hydrogen/oxygen  partial  pressures  within  the  corresponding 
channels  to  obtain  the  empirical  MEA  model.  While  such  di¬ 
rect  utilization  of  hydrogen/oxygen  partial  pressures  within 
the  channels  could  lead  to  error  in  thermodynamic  potential, 
such  error  can  readily  be  compensated  for  in  the  activation 
overvoltage  calculation.  A  detailed  comparison  will  be  pro¬ 
vided  in  Section  3.2. 

Generally,  fuel  cell  voltage  is  the  summation  of  three  ef¬ 
fects,  the  thermodynamic  potential  ^Nemst,  the  activation 
overvoltage  rizct  resultant  from  anode  overvoltage  and  cath¬ 
ode  overvoltage,  and  the  ohmic  overvoltage  77ohmic-  For  hy¬ 
drogen/oxygen  fuel  cell,  the  thermodynamic  potential  is  [19] 

^Nernst  =  1-229  -  (8.5  X  l(T4)(7body  -  298.15) 

+  (4.308  x  10-5)7body[ln(^H2)+0.51n(Po2)] 

(14) 

Here,  we  will  use  the  channel  hydrogen/oxygen  partial  pres¬ 
sure  instead  of  the  effective  partial  pressure  used  in  [9].  The 
activation  overvoltage  can  be  written  as 

7? act  —  §1  T  §2^body  T  §3^body  hl(/)  +  §4^body  hl(co2)  (1^) 

The  parameters  in  Eq.  (15)  can  be  obtained  from  experi¬ 
mental  data  using  a  linear  regression  technique.  While  var¬ 
ious  empirical  expressions  for  activation  overvoltage  exist, 
one  can  always  use  curve  fitting  to  obtain  a  proper  mathe¬ 
matical  model  [9,14,18].  The  ohmic  overvoltage  is  given  as 

-v.  _  ;  ointernal 

7? ohmic  —  lP  (lb) 

where  7?internal  is  the  lumped  resistance  caused  by  resistances 
in  electron  flow,  proton  flow  and  other  contact  resistances. 
The  empirical  equation  for  internal  resistance  uses  the  format 
as  follow  [9] 

Vernal  =  +  ^rbody  +  ^7/  (17) 


2.2.  Electrochemical  modeling 

In  previous  studies,  empirical  modeling  has  been  devel¬ 
oped  for  the  MEA,  which  is  built  upon  a  phenomenological 
relationship  between  the  fuel  cell  voltage  and  various  con¬ 
tributing  variables  [8,9].  This  approach  has  been  employed 
in  various  system-level  fuel  cell  modeling  [12-15].  Typi¬ 
cally,  the  cell  voltage  is  assumed  to  depend  upon  state  vari¬ 
ables  including  the  partial  pressure  of  hydrogen  and  oxygen, 
the  temperature  of  the  cell  body,  and  load  current,  etc. 

Vcell  =  /(Eh2,  Eo2,  Tbody,  U  Ot)  (13) 


where  the  coefficients  will  also  be  solved  using  a  linear 
regression  method.lt  should  be  noted  that  Eqs.  (15)  and  (17) 
apply  only  to  the  steady-state  case.  In  order  to  capture  the 
transient  dynamic  behavior,  in  this  study  we  incorporate  the 
effect  of  double  charged  layers  [12,16]  by  defining  the  acti¬ 
vation  resistance  as 

Ra  =  (18) 

l 

The  dynamic  activation  overvoltage  is  then  governed  by 

dnact  i  nact 
d  t 


C  RaC 


(19) 
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which  yields  the  total  cell  voltage  as 


3.1.  Benchmark  study  and  analysis 


^cell  —  -^Nernst  T  ^act  T  ^ohmic  (20) 

The  modeling  of  the  overall  fuel  cell  system  is  then  com¬ 
pleted. 


3.  Benchmark  simulation  and  experimental  correlation 

Outlined  in  Section  2  are  the  main  dynamic  equations  of  a 
fuel  cell  system.  In  this  section,  we  present  a  series  of  anal¬ 
ysis  results  to  verify  the  developed  model  and  to  study  the 
complicated  dynamic  interactions  between  various  compo¬ 
nents  and  effects.  The  mathematical  model  is  implemented 
using  MATLAB/SIMULINK,  as  shown  in  Fig.  2. 


We  first  compare  the  newly  developed  model  with  the 
benchmark  work  performed  by  Amphlett  et  al.  [12],  and 
demonstrate  the  improvement  in  the  proposed  work.  In  that 
study,  while  the  temperature  transient  behavior  of  stack/cell 
is  considered,  the  modeling  of  gas  flows  is  based  on  steady- 
state  equations.  In  fact,  the  model  developed  in  [12]  can  be 
considered  as  a  special  case  of  the  dynamic  model  we  have 
developed  in  this  research.  This  can  be  illustrated  by  the 
following  analysis. 

Consider  the  anode  channel  control  volume  continuity 
equations  given  as  (1),  (10a)  and  (10b).  If  we  assume  steady- 
state  flow,  we  have 


i 

"up,an(Ps  Pan)  ^down,an(Pan  Patm)  =  N  Mh2  (21) 

Zr 


Tstack 


Fig.  2.  SIMULINK  block  diagram  of  the  mathematical  model. 
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In  virtue  of  Eqs.  (4)  and  (5),  one  can  get 


^H2,an,in  —  ^up,an(/^s  Pan) 

(22a) 

^H2,an,out  =  ^down,an(/?an  —  Patm) 

(22b) 

i 

^H2,cons  =  N  Mu2 

Zr 

(22c) 

Substituting  Eqs.  (22a)-(22c)  into  (21),  we  can  obtain  ex¬ 
actly  the  analytical  flow  rate  relation  given  in  [10].  Similar 
assumption  can  be  applied  to  cathode  channel  control  vol¬ 
ume,  and  we  can  again  get  the  same  conclusion.  Therefore, 
the  flow  rate  equations  presented  in  [12]  are  the  steady-state 
special  case  of  our  model. 

We  then  compare  the  energy  equations  of  our  model  with 
those  given  in  [12].  In  that  study,  Amphlett  et  al.  defined  four 
energy  terms,  the  theoretical  energy  produced  by  electro¬ 
chemical  reaction  (gtheoX  the  electrical  energy  output  (qQ iec), 
the  heat  loss  rate  from  stack  surface  through  heat  convection 
(giossX  and  the  sensible  heat  terms  for  each  of  the  fuel  cell 
streams  (anode,  cathode,  and  water)  (gSens)-  Here,  we  again 
use  anode  channel  for  illustration.  Based  upon  the  notations 
used  in  this  research,  the  aforementioned  energy  terms  sat¬ 
isfy  the  following  relations 


^sens^ 


(V-n)d  A 


(23a) 


*/loss,a  —  ^conv,an(2body  Tan) 


(23b) 


d(mH2Ct;,H2^an) 
d  t 


—  ^sens^n  T  ^loss,an 


(23c) 


Table  1 


Parameters  used  in  the  comparison  simulation  [12] 


Parameter 

Value 

Parameter 

Value 

A^H2,an,in  (mol/s) 

0.0078 

VrI2,an,out  (mol/s) 

0.004 

Tan, in  (°C) 

23.5 

Tan, out  (  C) 

25.3 

No2,an,in  (mol/s) 

0.004 

V)2,an,out  (mol/s) 

0.002 

Tea, in  (°C) 

23.5 

Tea, out  (  C) 

38.8 

(M) an  (W/K) 

2 

(M)c a  (W/K) 

10 

(M) h2o  (W/K) 

50 

(M)st aek  (W/K) 

17 

MC  (kJ/K) 

35 

If  we  assume  the  steady-state  operating  condition,  the  left- 
hand  side  of  Eq.  (23c)  vanishes,  which  yields 

^sens^n  =  <7loss,an  (23d) 

The  above  relation  is  exactly  the  same  as  that  shown  in  [12]. 
Similar  results  can  be  obtained  for  the  cathode  control  vol¬ 
ume  and  water  flow  stream.  Clearly,  the  energy  equations 
related  to  gases  and  water  streams  presented  in  [12]  are  also 
the  steady-state  special  case  of  our  model.  We  can  then  con¬ 
clude  that  the  model  developed  in  [12]  is  generally  a  special 
case  of  the  present  work.  The  above  analytical  demonstration 
shows  that  the  newly  developed  model  will  yield  the  same 
simulation  result  under  certain  steady-state  assumptions. 

One  thing  worth  noting  is  that  in  dealing  with  the  sen¬ 
sible  heat  term,  Amphlett  et  al.  did  not  take  into  consid¬ 
eration  the  hydrogen  diffusion  through  the  electrode  [12]. 
Such  effect  is  included  in  the  current  model  and  can  be 
illustrated  by  using  the  data  provided  in  [12]  (shown  in 
Table  1).  The  load  current  changes  from  0  to  20  A,  and 
as  a  result  the  output  voltage  changes  from  38  to  28.6  V 


Fig.  3.  Effect  of  diffusion  on  transient  stack  temperature. 
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[12].  We  study  the  stack  temperature  history,  and  the  sim¬ 
ulation  result  is  shown  in  Fig.  3.  Clearly,  the  gas  diffusion 
will  affect  the  stack  temperature.  One  can  further  envision 
that  for  a  larger  current  situation,  this  effect  will  become 
even  more  significant,  as  the  consumed  gases  will  further 
increase. 


3.2.  Modification  of  empirical  MEA  model  development 

Amphlett  et  al.  [11]  developed  an  empirical  MEA  model 
predicting  the  fuel  cell  output  voltage,  activation  overvoltage 
and  internal  resistance  based  on  experimental  results.  Such 
approach  has  been  adopted  in  many  system-level  PEMFC 


Fig.  5.  Output  voltage. 
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Table  2 


Empirical  MEA  model  coefficients 


New  model 

Amphlett  model  [11] 

$1 

-0.97252 

-0.9514 

£2 

0.0034442 

£2 

0.00312 

£3 

9.3216E— 005 

£3 

7.4E-005 

£4 

-0.00019141 

£4 

-0.000187 

£5 

0.016046 

£5 

0.016046 

^6 

— 3.4715E— 005 

^6 

— 3.4715E— 005 

^7 

7.9565E— 005 

^7 

7.9565E— 005 

modeling  studies.  In  that  approach,  the  effective  surface  par¬ 
tial  pressure  must  be  calculated  based  on  the  estimation  of 
the  pressure  diffusivity  product  and  the  thickness  of  the  gas 
diffuser  portion  of  the  electrode.  The  accurate  estimation 
of  these  two  parameters,  however,  is  not  straightforward.  In 
this  research  we  propose  to  simplify  this  procedure.  When 
employing  Fick’s  law,  we  use  the  gas  partial  pressure  in  the 
fuel  cell  channels  to  directly  solve  for  the  seven  coefficients 
in  the  empirical  relations.  Using  the  data  provided  in  [11], 
the  model  coefficients  are  obtained  with  a  linear  regression 
technique  and  are  listed  in  Table  2.  The  coefficients  of  inter¬ 
nal  resistance  model  are  the  same  as  those  given  in  [11],  as 
the  partial  pressure  has  no  effect  on  the  internal  resistance. 
While  using  the  gas  partial  pressure  in  the  Nernst  equation 
could  lead  to  certain  errors,  such  errors  can  be  compensated 
for  in  calculating  the  activation  overvoltage.  Fig.  4  shows  the 
polarization  curve  comparison  under  the  same  conditions. 
The  result  obtained  through  the  simplified  approach  shows 
good  agreement  with  that  given  in  [11]. 

3.3.  Fuel  cell  dynamic  model  simulation 

In  this  section,  we  carry  out  a  detailed  parametric  analy¬ 
sis  to  study  the  transient  dynamic  effects  within  the  fuel  cell 
system.  The  parameters  used  in  this  part  of  study  are  listed 
in  Table  3.  For  illustration  purposes,  the  parameters  we  use 
are  collected  from  the  literature.  Due  to  the  unique  improve¬ 
ment  of  the  proposed  model,  certain  parameters  have  not 
yet  been  considered  by  other  studies  and  we  then  use  those 
that  could  be  representative  of  typical  PEMFC  systems.  As 
some  parameters  such  as  lumped  channel  volumes  are  small, 
stiff  differential  equation  solver  ODE23S  is  employed  in  the 
simulation.  The  time  length  for  the  simulation  is  36,000  s. 

Fig.  5  shows  the  stack  output  voltage  change  correspond¬ 
ing  to  load  current  change.  At  12,000  s,  the  current  rises 
from  10  to  50  A,  and  the  voltage  drops  from  29  to  18  V.  At 
25,000  s,  the  current  goes  down  from  50  to  35  A,  and  the 
voltage  goes  up  from  18  to  22  V.  One  can  see  that  within 
each  stage,  while  the  respective  current  value  remains  con¬ 
stant,  the  voltage  value  keeps  changing.  The  varying  volt¬ 
age  is  the  result  of  combined  effects  due  to  the  capacitance, 
time  varying  cell  body  temperature,  and  partial  pressures  of 
hydrogen  and  oxygen. 

Fig.  6  shows  the  temperature  histories  of  the  stack  body, 
anode  channel  and  cathode  channel.  The  temperatures  of 


Table  3 


Parameters  used  in  system  simulation 


Parameter 

Value 

Reference 

£Up,an  (kg/(s  atm)) 

3.6  x  nr6 

a 

k down ,  an  (kg/ (  S  atm)  ) 

2.2  X  10“4 

a 

Ps,an  (atm) 

2.4 

[12] 

Patm  (atm) 

1 

Well  known 

Un, in  (K) 

296.5 

[12] 

Tan, out  (K) 

296.5 

[12] 

cv,u2  (J/(kgK)) 

10124.71 

[19] 

Cp,u2  (J/(kgK)) 

14209 

[19] 

Mh2  (kg/mol) 

2.016 

[19] 

kup,ca  (kg/(satm)) 

3.6  x  nr4 

a 

kdown,ca  (kg/(satm)) 

2.2  x  nr3 

a 

Ps,ca  (atm) 

2.4 

[12] 

N 

35 

[12] 

Tea, in  (K) 

296.5 

[12] 

Tea, out  (K) 

296.5 

[12] 

c,,o2  (J/(kg  K)) 

662.2 

[19] 

cp,o2  (J/(kgK)) 

922 

[19] 

Mo2  (kg/mol) 

32 

[19] 

c„,n2  (J/(kgK)) 

745.2 

[19] 

Cp, n2  (J/(kgK)) 

1042 

[19] 

kconv,amb  (W/K) 

17 

[12] 

kCOnv,an  (W/K) 

2 

[12] 

kCOnv,ca  (W/K) 

10 

[12] 

^7bodyC/?,body  (J/K) 

35000 

[12] 

F  (As/mol) 

9.648456  x  104 

[19] 

ISHr  t  (I) 

1.196  x  10s 

Well  known 

Volan  (m3) 

0.005 

[13] 

Volca  (m3) 

0.01 

[13] 

$1 

-0.9514 

[11] 

^2 

0.00312 

[11] 

-0.000187 

[11] 

7.4  x  nr5 

[11] 

^5 

0.01605 

[11] 

^6 

-3.5  X  lO"5 

[11] 

^7 

8  x  lO"5 

[11] 

C  (F) 

2 

[18] 

Troom  (K) 

296.5 

[12] 

a  These  coefficients  are  so  chosen  that  the  maximum  consumed  fuel/gas 
is  satisfied. 


these  three  control  volumes  exhibit  similar  trends.  During 
the  first  12,000  s,  the  temperatures  slowly  increase.  The  rea¬ 
son  is  that  the  heat  generated  by  the  electrochemical  reac¬ 
tion  is  more  than  the  heat  transferred  to  the  surrounding 
atmosphere  through  convection.  During  the  second  stage 
(12,000-25,000  s),  the  output  power  of  the  stack  increases 
due  to  the  load  current  increase,  and  correspondingly  the 
heat  generated  by  the  electrochemical  reaction  increases  and 
overcomes  the  heat  loss  due  to  convection.  Therefore,  the 
temperatures  of  stack/anode/cathode  all  increase.  During  the 
last  stage  (25,000-36,000  s),  the  temperatures  drop,  as  the 
output  power  and  hence  the  heat  generated  is  reduced.  It  is 
worth  mentioning  that  the  temperatures  of  the  three  control 
volumes  differ  significantly.  This  clearly  demonstrates  the 
advantage  of  dividing  the  fuel  cell  system  into  three  control 
volumes. 

Fig.  7  shows  the  pressure  variation  within  the  two  chan¬ 
nels.  When  the  load  current  is  high  (50  A),  the  fuel  consumed 
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by  the  electrochemical  reaction  increases.  Correspondingly, 
the  gas  pressures  go  down.  On  the  contrary,  when  the  load 
current  is  low  (10,  35  A),  the  pressures  within  the  channels 
become  relatively  high. 

Figs.  8-10  show  the  change  of  mass  of  gases  (hydrogen, 
oxygen  and  nitrogen)  within  the  two  channels  due  to  the 
load  current  change.  During  the  first  stage  and  last  stage,  the 
load  current  is  lower,  which  leads  to  lower  fuel/oxygen  con¬ 
sumption  for  the  electrochemical  reaction  and  higher  mass 


of  gases  within  the  channels.  The  opposite  can  be  observed 
for  the  second  stage.  It  is  shown  that  while  hydrogen  and 
oxygen  decrease  due  to  load  current  increase,  the  nitrogen 
mass  goes  up.  This  is  because  when  the  oxygen  consump¬ 
tion  increases,  the  pressure  within  the  cathode  channel  goes 
down  and  more  air  will  enter  from  the  inlet,  which  causes 
the  increase  of  nitrogen  mass  as  nitrogen  is  not  consumed  by 
the  electrochemical  reaction.  The  inlet  flow  rates  are  shown 
in  Figs.  1 1-13. 
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Fig.  7.  Pressure  changes  within  the  anode  and  cathode  channels. 
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x  10  ~8 


Fig.  8.  Hydrogen  change  within  anode  channel. 


The  simulation  results  presented  in  this  subsection  clearly 
illustrates  the  complicated  dynamic  interactions  between 
various  components  and  effects  (temperature,  gas  flow  and 
capacitance,  etc.),  and  demonstrates  the  necessity  of  devel¬ 
oping  a  transient  model  using  the  control  volume  approach. 


3.4.  Experimental  correlation  and  verification 

To  further  validate  the  modeling  strategy  and  provide 
more  insight,  a  prototype  fuel  cell  was  designed  for  experi¬ 
mental  test.  The  single  cell  consists  of  a  Nation  membrane 


x  10 


Fig.  9.  Oxygen  change  within  cathode  channel. 
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x  10  "8 


Fig.  10.  Nitrogen  change  within  cathode  channel. 


sandwiched  between  two  carbon  fiber  paper  electrodes.  The 
surface  area  of  the  cell  is  6.25  cm2.  In  the  experiment,  pure 
hydrogen  and  oxygen  are  supplied  to  the  anode  and  cathode 
channels.  For  experimental  purpose,  a  temperature  control 
unit  is  employed  to  introduce  designated  cell  body  temper¬ 
ature.  A  series  of  system  input/output  data  (current,  temper¬ 
ature,  output  voltage,  ohmic  overvoltage,  inlet  volume  flow 


rates  of  anode  and  cathode  channel)  are  recorded  by  the  on¬ 
line  measurement  system.  Three  sets  of  experimental  results 
are  obtained. 

For  each  set  of  experiment,  we  adjust  the  load  current. 
Typically,  the  fuel  cell  operates  for  300  s  under  one  specified 
load  current  value,  and  the  time  interval  between  consecutive 
measurements  is  selected  as  15  s.  Therefore,  20  data  points 
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Fig.  11.  Inlet  flow  rate  of  anode  channel. 
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Fig.  12.  Inlet  oxygen  flow  rate  of  cathode  channel. 


can  be  obtained  under  one  specific  load  current  value.  We 
measure  the  corresponding  voltage  output  and  other  depen¬ 
dent  variables  such  as  inlet  flow  rate  and  ohmic  overvoltage, 
and  the  results  are  listed  in  Table  4. 

The  volume  flow  rates  listed  in  Table  4  can  be  converted 
into  mass  flow  rates  using  the  following  equations  [19] 

p  =  pRT  (24) 


m  ~  pV  (25) 

The  mass  flow  rate  coefficient  given  in  Eq.  (4)  is  a  key 
parameter  in  exploring  the  dynamic  behavior  of  the  gases 
flowing  through  the  channels.  The  upstream  mass  flow 
rate  is  estimated  using  the  traditional  nozzle  equation 
[13] 
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Fig.  13.  Inlet  nitrogen  flow  rate  of  cathode  channel. 
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Table  4 


Experimental  data  obtained  from  single  cell  test 


T  (K) 

Current  (A) 

Total  output  voltage  (V) 

Ohmic  overvoltage  (V) 

H2  inlet  flow  rate  (l/min) 

O2  inlet  flow  rate  (l/min) 

363 

0.50003 

0.77271 

0.00971 

0.20523 

0.20188 

363 

1.00022 

0.71458 

0.03375 

0.21075 

0.20431 

363 

1.50024 

0.6672 

0.05078 

0.21037 

0.20704 

363 

1.99961 

0.62471 

0.06748 

0.22179 

0.2097 

363 

2.49986 

0.58464 

0.08467 

0.22725 

0.21266 

363 

2.99999 

0.54541 

0.10229 

0.23278 

0.2155 

363 

3.50002 

0.50558 

0.12056 

0.23826 

0.21813 

363 

3.9999 

0.46499 

0.139 

0.24392 

0.2207 

363 

4.49976 

0.42315 

0.15851 

0.24924 

0.22368 

363 

5.00031 

0.37885 

0.17983 

0.25492 

0.22638 

363 

5.50012 

0.33093 

0.20306 

0.26036 

0.22891 

353 

0.50001 

0.81165 

0.00492 

0.20527 

0.20155 

353 

1.00003 

0.7755 

0.01429 

0.21086 

0.20411 

353 

1.50007 

0.75073 

0.02147 

0.2163 

0.20705 

353 

2.00029 

0.7283 

0.02955 

0.22178 

0.20989 

353 

2.50014 

0.70688 

0.03761 

0.22732 

0.21246 

353 

2.99981 

0.68706 

0.04482 

0.23286 

0.21518 

353 

3.50002 

0.67186 

0.04766 

0.23837 

0.21799 

353 

4.00011 

0.65476 

0.05586 

0.24384 

0.22091 

353 

4.49992 

0.63519 

0.06545 

0.24931 

0.22353 

353 

5.00035 

0.61905 

0.07251 

0.25478 

0.22647 

353 

5.50008 

0.60365 

0.07991 

0.26036 

0.22882 

km.up  =  6.3  X  10  12kg/(sPa), 
kca.up  =  7.5979  X  10“ 11  kg/(sPa) 

Clearly,  the  outlet  flow  rate  coefficient  depends  upon  the 
pressure  difference  and  temperature.  Here,  we  treat  the  tem¬ 
perature  as  constant,  as  the  temperature  is  indeed  kept  as 
constant  in  the  experiment.  An  nth-order  polynomial  is  used 
to  interpolate  the  outlet  flow  coefficient 


n 

kdown  =  y>(A/y-''  (26) 

i= 1 

The  coefficients  in  the  above  equation  are  obtained  using  the 
experimental  results  as  shown  in  Table  5.  After  we  obtain 
the  mass  flow  rate  coefficients,  the  pressures  of  hydrogen 
and  oxygen  in  the  channels  can  be  solved  using  Eq.  (1). 
Using  the  simplified  method  outlined  in  Section  3.2,  the 


Fig.  14.  Load  current  of  experiment. 
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Fig.  15.  Comparison  of  experimental  and  simulated  output  voltage. 


Table  5 

Outlet  mass  flow  rate  coefficients 


Anode  channel 

Cathode  channel 

ci  =  4.7433E— 030 

a  =  — 2.084E— 034 

c2  =  — 1.08 12E— 024 

c2  =  9.1435E— 029 

c3  =  9.7248E— 020 

c3  =  —  1.6348E— 023 

c4  =  — 4.3164E— 015 

c4  =  1.5164E— 018 

c5  =  9.0445E— 011 

c5  =  — 7.6298E— 014 
c6  =  1.886E— 009 

Table  6 

Coefficients  associated  with  MEA  model 


Coefficients  in  Eq.  (15) 

Coefficients  in  Eq.  (17) 

£i  =0.581 

£5  =  -0.133 

£2  =  0.02536 

|6  =4.103  x  10“4 

|3  =  1.044  x  10“4 

|7  =  7-9584  x  10“4 

|4  =  -1.0028  x  10“3 
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Fig.  16.  Comparison  of  experimental  and  simulated  anode  inlet  flow  rate. 
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Fig.  17.  Comparison  of  experimental  and  simulated  cathode  inlet  flow  rate. 


coefficients  related  to  activation  overvoltage  and  internal  re¬ 
sistances  can  be  solved  using  a  linear  regression  technique, 
which  are  shown  in  Table  6.  Two  sets  of  experimental  data 
are  used  for  identifying  the  key  parameters  used  in  the  dy¬ 
namic  model,  and  the  other  set  is  used  to  compare  between 
the  simulation  and  experimental  investigation.  The  load  cur¬ 
rent  history  is  shown  in  Fig.  14.  To  facilitate  the  simulation, 
the  experiment  data  are  interpolated  into  continuous  current 
signal  with  a  third-order  spline  function.  Here  the  current  is 
used  as  the  input  signal  to  the  dynamic  model,  and  the  output 
voltage  and  inlet  flow  rates  of  anode  and  cathode  channels 
are  simulated. 

Fig.  15  shows  a  comparison  between  experimental  and 
simulated  output  voltage.  In  general,  the  simulation  results 
exhibit  excellent  agreement  with  actual  cell  performance.  At 
very  large  current  values,  one  can  observe  certain  deviation 
between  the  predicted  and  experimental  voltage  values.  One 
reason  is  that  a  very  large  current  will  cause  a  large  num¬ 
ber  of  water  molecules  be  dragged  by  the  hydrogen  proton 
from  anode  side  to  cathode  side,  which  could  lead  to  the 
flooding  of  the  cathode  catalyst  and  membrane  drying  at 
anode  side.  Both  effects  could  actually  deteriorate  the  per¬ 
formance  of  the  cell  [1-7].  Another  possible  reason  is  the 
cell  temperature  effect.  Only  two  cell  temperature  values  are 
used  in  the  experiment,  80  and  90  °C.  These  two  tempera¬ 
ture  values  might  not  be  enough  to  cover  the  entire  operating 
range. 

Figs.  16  and  17  show  the  comparison  between  the  exper¬ 
imental  and  simulated  inlet  flow  rates  at  anode  and  cathode 
channels.  Again,  one  can  see  excellent  agreement  between 
model  prediction  and  experimental  result. 


4.  Conclusion 

In  this  study,  we  have  developed  a  system-level  dynamic 
model  of  PEMFC  that  incorporates  the  complicated  tem¬ 
perature,  gas  flow  through  the  channel  and  the  capacitance 
formed  by  double  charged  layer  of  MEA  under  operating 
conditions.  To  clearly  quantify  the  dynamic  interactions,  we 
divide  the  fuel  cell  system  into  three  control  volumes,  the  an¬ 
ode  channel,  the  cathode  channel,  and  the  fuel  cell  body,  and 
develop  respective  lumped-parameter  dynamic  models.  The 
system-level  model  is  simulated  using  SIMULINK.  Exten¬ 
sive  numerical  studies  are  carried  out  and  it  has  been  demon¬ 
strated  that  the  developed  model  can  capture  the  transient 
dynamic  behavior  of  the  fuel  cell  system.  We  have  also  per¬ 
formed  an  experimental  study  and  have  illustrated  excellent 
agreement  between  the  predicted  and  experimental  results. 
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